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Features of the computer program include the following: 
use of many subroutines for clarity and for easy modification, 
and use of program structure for easy integration into an 
optimization code. 

Sample lens designs and lens performance are presented. 
The computer code provides the tools for lens design. Consider- 
able additional work is needed for selecting the best lens 
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1. INTRODUCTION 


A long-term conflict exists between optics and aero- 
dynamics for missile nose design. Optics requires a blunt 
shaped nose similar to that used on the Sidewinder missile. 
Favorable aerodynamics for supersonic missiles requires a sharp 
nose which permits an attached shock wave. For rocket 
powered missiles or projectiles, the sensor and associated 
optics occupies the complete nose area. The penalty of a blunt 
nose is high drag. For ramjet propelled missiles or pro- 
jectiles, the penalty of a blunt nose is poor pressure recovery 
by the inlet; this statement assumes a nose inlet is used. 

In times past, the lens designer had available three design 
parameters; a satisfactory lens in the general shape of a cone 
has not been designed. A lens can be designed using homo- 
geneous refractive index which focuses rays; however, off-axis 
rays have very large aberrations. The thesis by Terrell [1] 
discusses the design of a conical lens. 

A fourth design parameter has become available which is 
Spatially variable refractive index. The associated technology 
is known as GRIN; see references [2-13]. GRIN is short notation 
for Gradient Refractive Index. The thesis by Frazier [14] 
briefly introduced the use of GRIN in lens design; likewise 
the thesis of Terrell [1] has a brief discussion. The thesis 
of Carr [15] will have more extensive results. 

The development of materials which have variable refractive 
index has been motivated by optical fiber technology. Optical 


fibers can propagate a single mode with low losses if an ap- 


propriate spatial variation of refractive index exists. Based 
on the information in reference [2] , the maximum variation 
of GRIN is limited to a few percent. The report contains 
twelve sections. In section 3 the problem of the homogeneous 
case is solved (program LENS) while the GRIN lens design is 
formulated in section 4 (program GRIN). The flow charts to 
the computer programs are given in section 5 and 6, and a 
user's guide to each of the programs is illustrated in section 
7. The programs are listed in sections 9 and 10. 
2. BRIEF OVERVIEW OF PROBLEM FORMULATION 

Define the axial distance as x and radial distance as 
r ina cylindrical coordinate system. The lens is to be 
axisymmetric so that angle is not a variable. 

In general the shape of the front surface, is a function 
Gf “x. 
(r) (229 
where subscript f£ denotes front surface. Likewise the rear 
surface is a function of r 

% = Ete) (22-2) 

where subscript r indicates the rear surface. 

As a Simplification, either the front or rear surface is 
a circular cone. Terrell [1] treated the case of a conical 
front surface; Carr [15] is investigating the case of a conical 
rear surface. This report develops a computer code which 
calculates either case. A single consistent computer program 


is used for both cases; this fact is important for application 


of optimization codes. 
In general, the GRIN can be specified as 
n = F(x,r) (233) 
The function F involves both spatial variables. However, a 
less complex GRIN is used which is spherically symmetrical 
refractive index. In this case 
n = n(x ,r ) (2.4) 
where r' is radial distance from the center of symmetry for 
the sphere. The location of the center of the spherical GRIN 
is on the lens axis at axial position Xue 
One very important computational advantage of spherical 
GRIN is that the rays within the lens are confined to a single 
plane. In the case of skew rays, the plane of the rays is 
not a meridianal plane of the lens. 
The equation for a ray trajectory in a spherically 


symmetric GRIN geometry is 


8, =o) + =| dx! (2.5) 
= rn (r‘)-E 


where E is a constant for a ray. The value of E is obtained 
from 

E = eng (rg) rgsin¥, (2.6) 
Define tr! as a unit vector in the direction of a line extending 
from x, to point P in Figure 1. Define R as a unit vector of 
a ray which has been refracted by the surface at point P. The 


angle Yo is obtained from 


(C9) 


cos Yy = rick (247) 


CENTER OF SYMMETRY OF GRIN 


Figure 1. Geometry for Ray Tracing in a Spherically Symmetrical 
GRIN medium. The ray vector, R , and the radius 
vector, r' , are illustrated. 


In equation (2.6), the € is a Sign function having a value of 


+ 1.0 or - 1.0 depending on the value of ¥ within the lens. 


ale (25 
(= 226 Y > n/2 
Equation (2.5) can be integrated for certain functions n(r'). 
One function is 
n(r') = [A + B(r'/r,)? | 7 (2: 


where A and B are constants. fy, is a value associated with a 
specific lens design; a convenient value for sa is that value 
for the ray defining the lens aperture. 

Combining equations (2.5) and (2.9) and completing the 


integration yields 


ae 
“af ee 4 2B“ /r',"-A 
6j(e') = -& | sin” ) —y»—s— sin 
| (a*+4BE“/r,*)* |(a*+4BE“/r)*) 
(2 


Equation (2.10) is an explicit relation for o. as a function 
of xr". ru: T, are the radii to the rear surface and to the 
front surface respectively. An explicit relation for r' as a 


function of Be can be derived; the result is 


v2 |B | 


z! (8) 
(2 ,an2;2 2 1 Sac uke 
A+ A°+4BE [To * sin ¥ 26 +sin 


Ol 


8) 


9) 


« LO) 


For details look at section 4.1. For certain phases of the 
calculation equation (2.11) is needed. 

Thinking ahead to numerical optimization, the function 
of equation (2.9) was chosen to avoid numerical integration 
of equation (2.5). Equation (2.5) must be evaluated several 
times for each ray. A lens design may involve 100 to 1000 
rays. Numerical optimization of the lens requires calculation 
of many different lens shapes. Hence, evaluation of equation 
(2.5) may occur several thousand times in the process of lens 
optimization. The GRIN-related parameters to be determined 
in optimization are A, B and xO- 

3. EQUATIONS FOR THE HOMOGENEQUS CASE 

3.) Meridian Plane 

The homogeneous case is included here for several reasons. 
First, the homogeneous case is a simple example in which the 
intricacies of lens design are apparent in considerable clarity. 
Iteration is not required for design of a lens with a homo- 
geneous refractive index, whereas for GRIN, iteration is 
necessary. Second, the homogeneous case provides a meaningful 
check of the GRIN program. To conduct the check, the value of 
B is zero in equation (2.9). Obviously the answers from GRIN 
program with B equal to zero must agree with.answers from the 
homogeneous case. 


The book by Kingslake [16] has been helpful in the program 


development. 
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3.1.1 The Equations 
Refer to Figure 3 for definition of the symbols. 
n(J) + I,' (J) = 8(I) (a5) 
sin(I,[J]) = N35 ° sin(I,'[J]) (352) 


The notation N35 is used for the ratio n3/n>- 


I,‘ (9) = (m/2-a) - (n(J) + T,(5)) (3).8) 
sin(I,[(J]) mig y * sin(I,'[J]) (344) 
a + I, (J) + U = 1/2 (355) 


3.1.2 Solution 
The angles a, 8(J) are inputs and are known. 
From eq. (3.5): 
I, (J) =71/2-a -U (3.5 5)) 
From eq. (3.4): 
I,'(J) = sin™*In,, + sin(1,[31)] (3.4.1) 


Substitute eq. (3.5) to eq. (3.3): 


Weta) ot T(J) = T, (9) = I,‘ (9) + U (323o)) 
But from eq. (3.2): 
Io'(d) = sin” *[n3, > sin(I,[J])] (36202) 
Promeq. (3.3.1): 
Substitute eq. (3.3.2) into eq. (3.2.1): 
Gy Seine ips > Simla = ty" (ay = nial UP] |= 482202) 
Substitute eq. (3.2.2) into eq. (3.1): 
Coes) 


Equation, (3.1.1) .-is “solved for 1(o):: 


n Sin (DZ, G) < Be I py —anitd)) +0), = sani sr = nj 
23 1 1 (3.5422) 


Use the formula: 

sin(A + B) = sinA cosB + cosA sinB 
And substitute: 

x = sin(n(J)) ; Yy _ 52 = cos(n(J)) 


Be 
Equation (3.1.2) will turn intozthe form: 
1 
x = / (33423) 
ep tene 
Where: 
n cos(I,(J) - I,'(J) + U] - cos[B(J) ] 
Di a ah (sa) 
N53 ° sin[I, (J) Ly oes U)] = sinisitoia 


Therefore: 


n(J) = sin ~[x] (31,5) 


Knowing n(J), I,' (J) can be solved from equation (3.1), and 
then, I,(J) can be solved from equation (3.2). 
3.1.3 Coordinates of the Front and of the Rear Surfaces 
3/1.3-1 The First Ray 
All the magnitudes will be normalized 


with respect to the radius R. Ris the aperture for the lens. 


Therefore: 
y,()) = R 


Ri = 2.0 R 
x, (1) 


(3.6) 
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The thickness of the lens is defined: 


T(1) = / (x3 ()) - cid + [y,() - yy (1) 1° (3.7) 


Therefore, the appropriate coordinates of the rear surface are: 


x. (1) = x, (1) + T(1) + cos[1T/2 -a - I,‘ (1) ] 
(oe8) 
¥Yo(1) = yy(1) - T(L) + sin[m/2 - a - T,'(1)] 


3eideSse2. The Other Rays: Coordinates of the Rear 


Surface 
The new point of the rear surface 
X5 (J) ,y5 (9) will be the intersection between the surface and 
the ray. 
(a) The equation of the surface: 
Y2(J) = m[{x,(J) - BB'] (3.9) 
m= tan[7/2 - n(J - 1)] = cot[n(J - 1)] (359.3) 
BB' = BF - B'F 


B'F = y5(9 = 1), * |cotin/2 = nt — 1) J+ cotl@@W —.1) ]] 


oS ytd = 1) + {tan[n(J - 1)] + cot[8(7 - 1)]} (32922) 

ee yo(J) = cot[n(J - 1)]{x,(J) - [BF - B'F]} (379...2) 
(b) The equation of the ray: 

Y5(J) = -tan[8(J)Jx,(J) + BF + tan [8 (J) ] (3510) 


ee y5(J) = tan[B(J)] - {=x (J) + BF} (3.9051) 


From equations (3.9.3) and (3.10.1), the coordinates of the 


next point at the rear surface are found; those are: 


s.(5) = BEltanie(g)) + cot(n(s = 1))} = cot{n(s-1)1_- B'P 
ae cotin(}] - 1)] + tan[8(u 


tan[8(J) ]{-x,(J) + BF} (311) 


kK 
i) 
ks 
tl 


at 


3.1.3.3 Coordinates of the Front Surface 
(a) Equation of the rear surface: 
¥; (9) = (tana) ° x, (J) (312) 


(b) Equation of the ray: 


y,(9) = m[x, (J) - X5 (J) ] ~ Y> (J) (3 e138) 
where: m= -tan[1t/2 -a - I,‘ (9) ] (3.3224) 
or: m= -cot[a + T,'(3)] (3523729 


From equations (3.12) and (3.13) the intercept point aes 


calculated: 
Veo) +", Cy  COeTG Te a. (a) 
2 Z ne 
x, (J) = 
tana + cot[a + I,'(J)] 
(3.4) 
¥;(9) = (tana) < x, (J) 
3.1.3.4 The Thickness at the New Point 
P(g) S"/ tah (d) Seas (a) 17 ey. Can eee (ays) 
Tg RD =i 74 Y2 : 
3.2 Skew Rays in the Homogeneous Case 
3.2.1 Coordinate Transformation 
The x,y coordinates will be shifted by the 
angle a. 
1 = ° —- % ° i 3} 6 dua 5 dh 
x x cos (a,,) 4 sin (a) ( ) 
ya sin (a,,) Pay ae cos (a,,) (3162) 
Zi =e (3.46. 3) 


3.2.2 Calculate the Intercept Point of the Front 


Surface 


Assuming that the (0,0) point is at point A 


v2 


(Refer to figure 2.), the following equations can be derived: 


(a) Ray: R72 =y? 


_ 
ee Ray = costa sy - (X + AB) ° oats 


Where R 
ray 


Zo! are known. 


(b) Cone: R = tan(a) 


cone 


Solving equations 


22 
Ng 
YO = 2 
cos (a 
\~ = Zo! 
Where: 
A,D, =. (tang) 
Yo’ 
CDS Sao tana 
p 
ee 


3.2.3 Direction Cosines of the Normal 


2 2 


faye 2 2 (rane) ~ (2 oe ee) 0 


of = Saeand Oe + AB) 


is the radius of the ray and the values for Yo'r 


for x results: 


6 I hy 


(ied Wasen) 


(3.19) 


Gs ro 


(35.1972) 


(33 19.53) 


(3s20.1) 


(S202) 


(3420.3) 


Spee) 


(Se215-1) 


Figure 4. Diagram showing coordinate transformation for 
Off-axis: rays. 


[az] = 2¢(tana)*(x + aB)? + y? + 27)7 


Substituting equation (3.21) can turm equation (2.21.4) 


the form: 


(x + AB) 


[xe _ 2tana 
cosa 


The direction cosines of the normal are therefore: 


_ (S£/6x) 


— = -sina 
[Az| 


oe 


= SS&/Sy) _ 2y 


= nea 2tana/cosa(*X + AB) 


Define: o = tan +(2z/y) 


y° ~ = (tana) 7 (x + AB)? 
rael + tan“c) = tan“a(x + AB) 


y = (coso) (tana) (x + AB) 
Therefore: 


ny 


Similiarly: 


(coso) (cosa) 


oy MO Oma 


? = (sino) (cosa) 
wy Weel 


3.268 Refraction at the Front Surface 


The skew ray R will be defined: 


R=R, - (e,) +R, * (ea) +R, + (eg) 


LS 


(Seelee) 


(321.4) 


into 


22a oO) 


(3.22) 


(34268) 


(352371) 


(302302) 


(Sadan) 


(3.24) 


Figure 5. Cross section of lens in a plane normal to the lens 
axis. Point P is the point where an off-axis 
ray intersects the front surface of the lens. 


Where R,, R,, R, are the direction cosines of the ray, 


and they have the value: 


= . = ae SF ° = 3. 
Ry cosa, : Ry s a : Ro 0 (3.26) 
The surface normal is defined by the following expression: 
as “a “~ “w 
i = ny ° (ey) + ny - (ey) ~ n, (e,) (3.27) 
Where the direction cosines of the normal surface were 


defined previously (equations 3.22, 3.23.3, 3.24). 


The refracted angle at the front surface is defined: 


1 R.N 
§ = cos a Sn Ie (3.28) 
in| . In| 
= ata ey 2 ae 
Where: |R| HR, - Ry + Ry = l 
CS e290) 
/ 2 2 2 
> / Pah _ 
li = /ny + ny + ny = 1 
Therefore: 
-1l 
= { } 3 
9 ccs LRy - ny, + Ry - ny + R, ny) (32.30) 
¢=T — 8 
= (3.31) 
a — 
¢' = sin [n)> Sind] 
3.2.5 Refraction Equations 
Referring to figure 6 the refraction at the 
rear surface with regard to the front surface will be as 
foliows [16, p. 147]: 
Bt > : 7 13 
noR' =n,R + [n,cos (¢ ) n,cos(¢) JN 
The direction cosines of the normal surface were defined in 


LT 


equations (3.22), (3.23.3), (3.24) and the direction 
cosines of the ray at the front surface were defined in 
equation (3.26). The direction cosines of the skew ray at 


the rear surface, will therefore be: 


N12 Ry + [{cos(¢') - Ny eos'(@)) |i .. ny 


LY = N12 Ry + [cos(¢') - N15 cos (¢)"] . ny (3.32) 


ls N15 Ry + [cos(¢o') - N45 cos(¢)] . ny 


3226 Intercept at the Rear Surface 
The equation of the cone of the rear 
surface is: 


y +2 = fax = 5b] (8533) 


where a is the slope and is defined from figure 7. 


R R - R 
a= tan (5) = re = a (32332) 
es | Vv J + i J 
b is the intercept at the X - axis and is calculated from 
equation (3.33.1): 
R oS ke Re el ok 
b=x, = Sti fst (2, 3am) 
J+ i J 
(XI, YI, ZI) is the point of intercept atethe rear surface. 
The equation of the rear surface is therefore: 
2 2 
YI" + ZE~ = (axt — b)) 50 (3. 3830 
where a,b were defined in equations (3.33.1), (3.33.2). 
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Figure 6. Diagram used to calculate the direction cosines 
of skew rays at the rear surface. This diagram is 
based on Kingslake [16]. 
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Figure 7. Geometry of a truncated cone which represents a 
small portion of the rear surface of the lens. 
The cone angle is es: 
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The direction cosines of the ray at the rear surface were 


defined previously (3.32). They are: 


Kt = 42 =_80 SE ao > Me Ze: 


(XO, YO, ZO) is the point of intercept at the front surface 


Pee ly el Ie ey Se t9. 3). 


/ 
DS al(et = x0) + 4vne= ¥o)- #<(2r =. 20)- 


Solving for D from equations (3.33.3) and (3.34) one gets: 


pe ek ea aess 
at a’ 
Where: 
BY 2 tye + iMPy* = actKt)- 
b' = (L')YO + (M')ZO = (K')a* (XO —- b/a) 
Bie yor ee =) a RO Soba 


Knowing D, the intercept point at the rear surface can be 
verified using equation (3.34): 

XI = K'D + XO 

YI = L'D + YO 


ZI = M'D + ZO 


ee Pe Equation of the Normal to the Rear Surface 


(3535) 


(32.36% 


(35507 


(oesbie 


The equation of the cone of the rear surface 


was defined previously (eq. 3.33): 


Site ye ee ea Gy ace ace SO 


1) 


2) 


3) 


-1) 
ae) 


~ 3) 


a,b areidefinedtin equations (2533 .4)5 (3.38 22)i< 


Vf = ey et Gy eet Se oe (3.3822) 


Similar to what was done at the front surface: 
== = -2a(ax = b) (3.38232) 


CS ie Z : 6f _ 
Yi); SZ 22 (3238550 


Hy 
aot 
it 
ml 
Da 
Ith 
N 
+ 
it 
+ 
=~ 
ob 
ro 


= 21 arae - aie ~ ee + al i 
oie =) td aed (33802) 


The direction cosines of the normal to the surface are as 


follows: 
_ (S:£/6x) -a 
n = —-—— = (3.39.2) 
Ae gs, ite (a eo 
ny = 5 ee (3339.2) 


(S52) 0h eae 


Z 


n. = 
es e 
ae 2 ‘oie 


Seces Refraction at the Rear Surface 


Similar to what was done at the front surface: 


K' (ey) + L' (ey) + M'(e,) (3.40.1) 


- 
u 


ny (ey) + ny(ey) + ng (e,) (3.40.2) 


> Saami. 


Where the direction cosines of the refracted ray (R') and 


> 
N 


the normal (N') were defined previously (equations 3.32, 3.39). 


The refracted angle at the rear surface: 


al Ri oN! 
pauses hi eee (3.41) 


os = Kin, + E'n. + Min (a\e4a),-1) 


eel . 
sin [n53 : sind, ] (3341.2) 


<>- 
t4 
i] 


3.2.9 Direction Cosines of Rear External Refracted 


Ray 


As done previously (equation 3.32): 


R= n53K' on [cos (o,') - n53c0s (¢,) ] ny 
ge N53L' - [cos($,") - n4,c0s (o>) ] - ny (3.302) 
a= N52! - [cos (¢,") ~ n,c0s (>) ] Ne 


3.2.10 Intersection with the Image Plane 


Xo = BF = DK ae ae 

Vv = " Vv 4 3 

Yom DL + YI (3543) 
= A Ww 

Zon DM + ZI 


Solving equation (3.43), the intercept with the image plane 


can be calculated: 


D, = (BF - XI)/K" 
vy. = (BET XE) ae 4 yt (3.43.1) 
Im K acer 
_ (BE - XI \ uu 
ba = (gE) mt + 
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4. FORMULATION OF EQUATIONS FOR GRIN LENS DESIGN 


4.1 Solving the Integral for or 


As mentioned in section 2: 


rls ni) )=. EY - 


a6 
rE dq. 
j63, 08s =| Sy kine eee i-2) 
oi 


Les ae are the radii to the front surface and to the rear 


surface, respectively. 


| pil Ser no (r-) -©p- sinv , (4.2) 
n=,/A + (i) (4.3) 
Substituting: 
V = & ' ; = -2r,*r (42a 


Equation (4.1) turns into the forms: 


V 
Bo 
AGO ee oe qv (4,62) 
e 0 2Xy Iz 
v AV + b- V se 
f£ ro. 
: 2 
Therefore: 2 (rp/r_) 
2(B/ 2D) WN = A 
sie oe 0 
Saladin aieameee Wir Ge 
(AW + ABEe7 oo) 2 
(r9/r¢) 
rc 2B J oe RE. a a 
_ smell - = eee £ 
a -4;/sin ie: uDEER OPTS FE Sk Ee (4.1.3) 
(A° + 4BE°/r,°) %  _ (AT ABRS a7) 


The appropriate expression for the radius will therefore be: 


7 V 2/E| 
eet) 5 


soos | was one) ee - A 
VA +/A + 4BE /Xo - Sin ¥26 + Sin 7 aa oe a 
JA” + 4BE /XQ 
(4.4) 
If the radius to the rear surface (xr) is calculated froma 
known radius to the front surface (r¢) the minus sign is used 
in front of the angle 94: This situation occurs with skew 
ray calculation. The plus sign is used when ie is calculated 
from a known rp: This situation occurs for lens design. 
4.2 Meridian Plane 
4.2.1 Solution for the First Ray 


(a) Define Point X, (1) ,¥, (1) 


As in the homogeneous case: 


R= 0 

¥, (1) = R (4.5) 
7 R 

X, (1) ~ tana 


(b) Define Angle to the Front Radius 


R 
tan(6, + 6.) = (4.6) 
0 e x 1) - Xa 
re(l) = 4) = — uit (4.7) 
sin (> (1) + 9, (1) 
Refering to figure 8, 99 is the angle to the rear radius, a, 
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is the angle between the rear and the front radii, and Xo 
is the location of the center of symmetry. 


(c) Choose es from Initial Thickness 


8 (1) T (1) /r,(1) 


99 (1) (8 9 (1) + 8 (1) ) - 601) 


(dad) Solve for Ve 


a + T, (1) +.u = 7/72 


I, (1) =W/2-a- u 


Knowing r,(1), the initial index of refraction (n,(r,)) 


can be calculated from eq. (4.3). Therefore: 


sin(I,(1)) = no(r_) /n, : sin(I,'(1)) 
ees me gin“ [n,/n.(r,) . sin(I, (1) 1] 


Refering to figure 8, Ve can be calculated. 
9, (1) + 84 (2) =a +t 6 
oe § = 6 (2) + 8, (1) - a 

6 + ¥,(1) - T,' (1) = 1/2 

v, (4) = at 2. tek 1G) 8S 


(e) Solve for E 


BS ess n,(X¢) s) 8 sinw 
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(4.8) 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


(42) 


*aATQeHSU ST q PpaPeARPSNTTT 9sed 3yF AOGT “*SuST NINO 
ayy Aq poeqz}tusueriq AeT e TOF Sa }eUTpIOOD pue saThue Huturhsep werhbetgq 


JDWJON 


Any 


DW JON 


°g omnbTtq 


™ 
ON 


+1 


WS 2a oo 
Where: e Sots (42258 


W> 1/22 -1 


(£) Solve for reid, 


Knowing the radius to the front surface 
(r¥) and the angle between the two radii (8.), the radius to 
the rear surface (r.) can be calculated using equation (4.4). 


Using equation (4.2), ve, can be calculated: 


ain 4 (4.2.2) 
r 


(g) Solve for 8 (1) 
The focal length, which is an input, is known: 
BF = F/R 


The coordinates at the rear surface are: 


Yj(1) = r-sin(9,(1)) | 
(4,13) 
Xx, (1) = Xo + r cos (8) (1)) 
Using equation (4.13), 8(1) can be calculated: 
ay Y, (1) 
B(1) = tan [seer | (4,14) 


(h) Solve for the Refractive Angles at the 
Front Surface 


As in the homogeneous case: 
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n(3) + In"(J) = (9) 


sin[I,(J)] = n3/n, (xr) . sin[I,'(J)] 


¥() - I, (J) + n(J) + 84 (od) = 


Solve for n(J) in a manner similar to the homogeneous case: 


n(g) = eso EG | 
Siler a 
siné' - n,/nj(r,) . sing (J) 
Where: 
6' =7T - 84 (9) - Vy (o) 


Knowing n(J), T,' (J) can be solved from eq. (4.15) and then, 
I, (J) can be solved from equation (4.16). 
4.2.2 Solution for Next Rays 
4.2.2.1 Coordinates at the Rear Surface 
Solve for the subsequent rays the 
same way as in the homogeneous case (equations 3.10, 3.10.1) 
4.2.2.2 Calculate Radius and Angle’ to the 


Rear Surface 


ae ae ee. 
xr (d) = / (Xo (I) = Sele 7 Y5 (J) 


ee: 
84 (J) =. tan [¥5 (J) /(X, (J) - Xa) ] 
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(4.15) 


(4.16) 


(4.17) 


(4.18) 


(4216.1) 


(4.18.2) 


(4.19) 


(4.20) 


4.2.2.3 Solve Front Surface Parameters 
Refer to Figure 8 for a definition 
of geometrical variables. 


(a) Equations 


n(J) + I4°(3) = B(d) (4.15) 
sin[I,(J)] = n3/n, (rx) A sin[I,' (J) ] (4.16) 
Yo Ad) F) Todd) ete nel gee tsar (4.17) 
) ae n,(r) By Sse at sin{y (J) ] (402.2) 
E. Sets no (r-) - rp - sin[).(J)] (4.2) 

Q,(J) + 6) (J) =a + 6 (4,21) 
é§ + ¥-(d) = T,' (J) = 1/2 (4.22) 
re = f(r, &,(9)) (4.4) 

sin[I,(J)] = no(r,)/ny, : sin[I,'(J)] (4.23) 
oe Ly (a) aS 3 7/2 (4,24) 


(b). Solution by Iteration 


An explicit solution is not 
possible. An iterative solution is necessary. 
(b1) Assume 8(9). 
(b2) Solve for coordinates 
of the front surface: 
- Calculate T, (3) from eq. (4.24) 
- Calculate § from eq. (4.21) 


- Calculate front surface coordinates and radius: 
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The 


Y, J) = tan(6_ (J) + 9, (I)) 5 (X4 (J) - Xa) C4225) 
Y, (9) = tance. X, (J) (4.26) 
result is: 
Xa . tan[6) (J) + 9(9) J 

X, (9) ~ tan[6, (3) ~ 6637] - tana (4.27) 

YX; (9) = tana . X, (J) (4.26) 

re = ¥4 (J) /sin[@, (J) ~ 8, (9) ] (4,28) 

Xa - tana 

= sin[@) (J) + Coe - tana . cos [8 (5) + 66 (dN) 
(4,526,.1) 


(b3) Knowing r,, I,'(J) can 


be solved using equation (4.23). 


for 


the 


- Calculate Ve from eq. (4.22) 
- Calculate E from eq. (4.2) 
(b4) Calculate second expression 
front radius (x) from eq. (4.4). 
(b5) The iteration routine: 
The iteration search 
correct o2 uses the Newton-Raphson routine. 
function: 


F = 0 (4.29) 


ree ye  e. 
is the expression from the formula (eq. 4.4) where r, is 
geometric expression (eq. 2.28). 


The derivative is found with respect to Pe 
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eles medic 
ff £ 
Flo = (4.30) 
a6.) a(@ 2) 
d(r,,) (zie -1 
SEE Bye : Bieelicesi=20. + sin (B5)] (4.30.2 
a —e ogee : 
e 2E 
Where: 
B. =/ar # ABleye Ne (4.30.2) 
al r e e 
z 2 


The lower signs are used when the front radius (x) is calculated 
from a known rear radius (x) (this case occurs at J > 1), while 


the upper signs are used when Bs is calculated from le (for J=l1). 


d(r,) “X, . tana. {cos [8 (J) + 6,5) ] + tana . sin[8) (J) + 6 (J 
até) Isin[6,(3) + 8,(J)] - tana . cos[@,(J) + a= (ae 
(4s, 30). 4) 


(c). Correct ‘Solution 
After knowing the correct value 
for afd), the procedure used in (b) can be followed resulting 
in the correct values for: 
T,(d), 7 X, (J), X09), re(J), T,'(J), ve), E. Using 
equation (4.2.2) Ue can be calculated. Equations (4.15), 
(4.16),(4.17) can be solved for n(J), I,'(d), 8(J) as shown 


previously for the J=1 case (equations 4.18, 4.15, 4.16). 


az 


4.3 Skew Rays 


4.3.1 Front Surface 
As in the homogeneous sace, the intercept 
point at the front surface (XO, YO, ZO) can be found (equations 
ee, 361932. 2el9~3) 3 
Also, the direction cosines of the normal to the front 
surface can be calculated (equations 3.22,3.23). The direction 
cosines of the ray can be found as well using equation (3.26). 
4.3.2 Front Radius 
Knowing the intercept point at the front 


surface, the front radius can be calculated: 


r_(J) = [(X0 - x)? + (vo)? + (20717 (4.31) 


Knowing r,(J), the initial index of refraction [n,(r_)] can 
be calculated using equation (4.3). 
4.3.3. Refraction at the Front Surface 
As for the homogeneous case (section 3.2.4), 
the refracted angle at the front surface can be calculated 
using equations (3.30), (3.31). 
4.3.4 Initial Direction Cosines for the Internal 
Refracted Ray 
The initial direction cosines of the internal 
refracted ray (K',L',M') are calculated the same way as for 
the homogeneous case (section 3.2.5, equation 3.32). The 


vector form of the initial direction of the ray is therefore: 


R= eK oe fi & eel Sout ce (44.32) 


4.3.5 Front Radius Direction Cosines 
The vector form of the front radius is 


defined as follows: 


Pe = fey = ey + rey ° ey + leg ey 
Where the direction cosines are: 
Tey (Xo - XQ) /¥_ (I) 
Fey = YO/r , (J) 
Yep = ZO/r - (J) 


(4.33) 


(4.33,1) 


4.3.6 Normal to Plane of the Curved Ray at (XO, ¥O, Zon 


The vector cross product of r. and R is normal 


£ 
to the plane defined by the two vectors. 


ey ey ey 
P.=r,.XR= 
NP Ee Ese t sy ee 
K' Gi M' 
NP; = NP cy : ey ~ NP ey ey + NP ey 5 en 
Where: 
= Taal ' 
NP 5y ley - M leg L 
= gone ’ 
NP py = Leg K lex M 
= LIN ape r] 
NP ep Tey L Ley K 
and 
rs 2 2 2 
NP, = NP py = NP ey = NP 29 


The normalized direction cosines of the normal to the plane 


of the curved ray at (XO,YO,Z0) are therefore: 


NP -y/NP - : NP .,/NP . ; NP -7/NP ; 
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(4.34) 


(4.34.15 


(4.34.2) 


(4.34.3) 


(4.34.4) 


4.3.7 Calculating the Angle Between the Ray and 


the Front Radius 


> 


ee [R | ) 
ase ° — IT ° | | e cos (Uv - 


Both the ray and the unit vector for the front surface radius 


are normalized, i.e. 


~ 2 
F,.| = ee yt HA) = 1 
Ia] = cet? tut? 4 my % 2 
Therefore: 
wv. = span ole R] 
if: if 
= coe e K' +r 
ce eed £Y 


4.3.8 Calculating constant E 


Knowing ., constant E can be 
using equation (4.2). 


4.3.9 Intercept Point of the Curved 


(4.35) 
(423564) 
T+ Leg - M' J 
(a e35¢2) 
calculated 


Surface (Point XI,YI,2I1) 
4.3.9.1 Equations 
As previously, 8. is 


the front radius and the rear radius. 


one 


x FX Ye Sale EY. rZ £Z 


The direction cosines of the front radius were 


the angle between 


= cos[@ (J) J (4°. 36) 


calculated in 


equation (4.33). Similiarly, the direction cosines of the 


rear radius are: 
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ba 
| 


(XI - X/ r(d) 


EX 
Coy = YI/r (J) (4. 326). 2) 
Eig ZI/xr (J) 


Therefore, equation (4.36) can have the form: 
(XI - Xa) Tey ~ (YI) rey - (ZI) Leg = r (do) : cos[6. (J) ] (4.536223 
The rear radius is defined as: 


(XI - i + (¥I)? + (21)7 = [x (3) 1" (4.37) 


The rear surface was defined previously. Consequently: 
S(ekt- seb). & er)- + ze S10 (3.33) 


a,b were defined by equations) 3.33.2 and 3.33.2. 
Also, from geometry: 


(XI - X) = rd) Z cos[6, (J) ] (4.38) 


Another equation arrives from the connection between the two 
radii through eq. (4.4). 
4.3.9.2 Solution 
a. Assume 9 (9). 
b. Calculate r (J) from the formula 
(equation 4.4). 
c. Calculate XI from equation 


(42 36%.2)-: 


SRS r (do) ; cos [8, (J) ] + X (4.36.3) 


@ 


d. Calculate radius from equations 
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(4.37) and (4.36.3): 


Radius = R = [ (yz)? e (zr) 7]? 


2 
r (dg) [1 - cos [89 (5) 1} 


= r (dg) - sin[@)(J) } (2237, 


e. Calculate parameters 


a and b from equations (3.33.1) and (3.33.2) as in the 
homogeneous case. 
4 
£. Substituting into equation 


(3553) 
AI = (R+b)/a = {fr (gd) . sin[6)(J)] +bi/a (4.38) 
where R is the radius defined in equation (4.37.1). 
g. Solving equations (4.36.3), 


(4.38) for 94 (9) one gets: 


5 ee vip'y? = (a') (c") 
AT a 


8, (5) = sin (4,39) 
Where: 
Aes ina a 4 
eos _ (4.40) 
B' = [xr (J)](b - a. X,) 
a) a p 
see (5) ee aa, Xa) , ae [r (J) ] 


h. Knowing 89 fd), it is pos- 
sible to calculate new value for 9 (9) and iterate through steps 


(a) - (g) until the difference in the values of 8) (J) is 


oe 


small enough. 
i. When an accurate value of 
84 (9) is received, XI is calculated from equation (4.36.3). 


YI and ZI are calculated using equations (4.36) and (4.37.1). 


Consequently: 
B" Jit) 2 pt CAS) (en) 
YI = ar te Mat oe se 
= (4.41) 
ZI = /fr (a) . sin[@,(o)1}° - (yr)? 
Where: 
eee: 2 2 
A" =1+ [rey] /(reg] 
B" = {r (J) ‘ cos [8 (J) ] - (XI - Xa) Ley} Pey/T ee 
rae ES) 1s 
Gap See leccipamanie=) (r= & oe : (3) . sanbeeeoia 
2 Ce sa 0 
ba £Z 
£Z 
(4,42) 
When Teg = 0, the solution for YI and ZI is different: 
= as ae + 
YI = uo) 7 cos[9.(J) ] (Sr Xa) Fey} /Tey (4.43) 


While ZI is calculated by equation (4.41). 

4.3.10 Direction Cosines of the Ray at the Rear 
Surface 
a. Similarly to what was done at the front 


surface (equation 4.35): 
; IR | : cos (Wy) (4.44) 
Both the ray and the rear radius are normalized; therefore: 


r el ee a a te lee vd Fe ee cos). (4.44.1) 
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b. The ray remains within the same plane 
while transfering from the front surface to the rear surface. 


NP - K_' + NP - L_' + NP - M.' = cos(t/2) = 0 (4.45) 
xr = r 


£2 


c. The direction cosines of the ray are 
normalized, i.e. 


12 12 2 
(Ce wy yes (4.46) 


ad. Solution: 


(4.47) 


w 
i 
ry 
Z 
rd 
fh 
K 
iH 
2 
' 
z 
rd 
Mh 
| 
fe 
H 
td 
Fh 
> 


Where: 


(4.48) 


Poy ) + 1 


2(NP cy) (NP ¢ ) 


Z 


AY = (at) + (Bt) (B)* - SEUSS) 
ree) SAB) rene OG) 
B" = oar 0:4 Bal (4.48) 
(B) 
o ee? 
Cree: 78 ) (=) - l 


In (4.47) the (+) sign is used when r < 0 while the (-) 


EY 
Sign is used when r = Or 

(a2) In the above case, when B = 0 (look for 
definition of B in equation (4.48)) the solution for direction 


cosines of the ray is different. 


a ee 
(L') == 
/ te") (L")) z 5 
(es) (L" ) yi Sa  = -(B") [(A') i =i) 
(Mt) 2 eo ee ee 4 NE a ere 
r 2(B") B! 
(K") = [-NP .., (L') - NP 5p (M') J/NP py 


(4.49) 


All the symbols used here were defined in equation (4.48) and 


the + sign is used as before. 


(M'_) = 0 (4.50) 
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ar 
/ 2 
(K! ) = B" n / (B") va (Am) (Cr) 
re ar a 
(4.50) 
(L') = — [cos (p_) - Ty: Ke 
ry = 
The (+) sign is used here. 
While in this case: 
beg ee 
A" = 1+ (===) 
Pry 
B" = cos (Ww ) rx ; 
(ry) (4551) 
Z 
; (eee) 
Crs = - l 
Ty, 
(d4) In the above case, when rey = 0 the 
solution is: 
(M re = 0 
cos (_) 
(K") = = 
~rX ; (452) 
’ ; / ’ 2 
(L' =i,/j/1- (K") 
Again, the (+) sign is used here. 
(a5) NP oy = 0, NP zy - Os NP <7 = 0, r7 = 0 


cee ee eee ae ee er ee es a a ee ee ee 


The solution is similar to that in case (d3). 


16 
B" J(B") 2 = GA) Cem) 
(K") = =r + ee (4:553) 
(Chee Be ae {[cos(W_) - xr Ke 
x sg re ap ee r 
eZ 


(A"), (B"), (C") are the same as defined in (4.51) except that 


08 is replaced everywhere by r The (+) sign is used here. 


ry ay Am 


(d6) In the above case, when E>. = 0 the 


solution is similar to that in case (d4): 


(L a) = 0 
c 
1 Cosh) 
g r x 
7 (4.54) 
' = / = ' 2 
(M » = tvJl (K =) 
Also here, the (+) sign is used. 
(d7) NP py = 0, NP cy #0, NP yy ~ 0, Ty = 0 
cos (y_) 
(L'.) =- est 
16; NP ey 
Try” eer TrzZ 
NP 
(M! ) = £Y Gee ) (4.555) 
1g NP cp 56 
ed ee eee 
(K =) Pino (L =) (M 2) 


The (+) sign is used 


(d8) 


(d9) 


As in case (dl), the 


the (-) sign is used 


here. 
NP 
, eV, 
As in (d7), but: (ry ip 
_ eZ 
Trx Ze 
cos (_) 
Trx 
/ 
re 1 \2 
i a (K 2 
PAM ies Se Nees 
/ 1+ (ese) 
ae £2 
os ae 
NP eo G 
As an “(d7),.but: rey 7 0 
w } w 2 we "Ww 
S(B" Dg AB )S = AA) (C™) 
yOu) a 2{A") 
is ae 
~ ~ NP (L' 
£Z 
NP 
1 £Y 
= —— [cos(v_) - (r 
rex 1g ry NP eg 
(+) sign is used when: ty 


. > z 
when: ry 2 0 
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< 0 while 


(4.56) 


5. FLOW CHART OF PROGRAM LENS 
(Homogeneous Case) 


Note: Subroutines SHAPET, SKEWT compute the shape of lens 
and the skew rays as described in the previous sections (3,4). 
However, subroutines SHAPEC, SKEWC compute the shape and the 
skew rays of a different design of the lens. In this design, 
the fixed cone is the rear surface instead of the front 
surface (compare figure 9 to figure 2.). The equations for 
this different design are not included in the current report. 
The equations are part of Capt. Carr's thesis [15] and will 


be described there. 
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Program LENS - Main 


DATTA 
YES 
| 
| (e) 
| | 
|  SHAPET SHAPEC 
| | | 
} | 
| 
| | 
| PRINT 
Change 
YES 
Gteeae) 
8 (J) 
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| 
| 


INTENS 


NO 


| ZOYO 
| 


RADI 


YES 


SKEWC 


| 


/ 


/ 
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/ SHAPE, 
SQURAY, 


ELLIPS 


(*) In program GRIN, which uses the same subroutine, here, 
data for XC, DIV is given. 
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SHAPET 


Compute 
T, (9) 


Compute 


1,‘ (9) 


Compute 


Z£,(d) 


| Compute 
| 5 Bea 8 


| 
Compute 


S| gta 


ops 


INIT 


Compute 


y,()) 


Compute 


x, (1) 


Compute 


Y, (1) 


Compute 


B (1) 


Gua 


POINTS 
| Compute 


x, (J) 


Compute | 


| 
{ 


| y., (J) 
an ae 
| Compute 


Compute 
| 


T(J) 


Compute 


nN (J) 
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SHAPEC 


| Compute 


1,1 (3) ,T, (5) | 


Compute 


| 
| 


erivertr | 
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Compute 


x, (J) r¥ (9) 


Compute 


£7 CS) ye 


1 69) 


Compute 


56 


Compute 


¥, (I) x, (9) 


Compute 


(*) 


* 
Compute | Compute 
AB ,BF AB 
Shift: Compute 
x1 (I) x, (9) Statna 


Compute 
y 


NO 


Compute 


Onp 


In program GRIN, which uses the same subroutine, here, 
also.xe 2s shafted: 53 


Compute 
Cen, Torbi, 
Boti 


5 nels 
Boti 


Compute 


Topi 


Compute 


¥6 (12) 


60 


Retum 


YES 


YES 


NO 


61 


62 


63 


Compute Compute 


R= Radius| Radius 


YES ICHECK = 40 


Compute 


Torb 


Cen, 
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Compute 


x0, yO, zO 


Compute 


Ni Ne No 


Change sign 


C6 


Compute 


Do 


Compute 


xi, yi, ZL 


Compute 


Ny 


Compute 


Orr ope ot 


YE 
: ICHECK=52 


no | 


68 


SKEWT 


Compute 


x0, yO, Zo 


Compute 


Koy, elt 


eS 
Compute 
D5 


ICHECK=51 


Compute 


O.: Op, ot 


7B 


Compute 


Ks ; yet P mM" 


Compute 


You (22) 2, 


12 


Compute 
xDiapt (Iz) 


xDiapt (Iz) =10 


yDiapt (Iz) yDiapt (Iz) =10 


Flag (Iz)=0 


Ntncty (Iz) =1. 


De fine 


r'(Iz “u 
yg {t ) ,25 (Iz) 


Count 
=Count + 1 


ICHECK=0 ICHECK=0 ICHECK=0 


Ntncty (Iz) =0. Ntncty (Iz) =0. tnety (Iz)=0. 


Flag(Iz)=1 Flag(Iz)=2 Flag (Iz) =3 


: 
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Flag (Iz) 


=Flag(Iz.) 
< 


Count=Count 
+] 


z}(Iz)= 
-z'i (Iz 


to 


xdiapt (Iz) 
=xdiapt (Iz.) 


ydiapt(Iz) 


=ydiapt(Iz.,) 


Op1 (Iz) 
=Op1 (Iz.) 


Ntncty (Iz) 
=Ntncty (Iz,) 


=< 
@> 


Flag(Iz.) 


So) 


62 


77 


Rays=Ia-1l 


78 


PRINT 


Coord=Coord 
+1 

El ()=E1 +El 
z ZZ 


YES 


a2 


El_()=El 
Zz ZZ 
-Ell 


Compute 


Coor()=Coord 


Coord =Coord 
+1 
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NO és 


YES 


Numbt=Numb-1 


a YES 


lag(Iz)> 


Sum = 
Ss 


Yom (t2) tSumy 


ai = .10 
angle p 


mhick=T (J) 


Compute 


R tead” 


Compute 


Sdrad (Iz) 


Compute 


Fractn (Iz,) 
i 


PRIN 


R 
oc(iz,) 


=Roc (Iz_) 
+0.005~ 
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85 


PLOTT 


PRINT (83) 


PRINT (84) 


PRINT (86) 
(87) 


eet YES 
(89) 


NO 


Coort= 
Coor (Numb) 


PRINT (97) 
(98) 


6. FLOW CHART OF PROGRAM GRIN 


Note: Program GRIN is identical to program LENS (Section 5) 
except for subroutines: 

SHAPET, SKEWT, SHAPEC, SKEWC. 
The flow chart for the first two is listed below. Subroutines 
SHAPEC, SKEWC are of a different design (fixed cone o.,, bale lrg 
the rear surface). They are part of Captain Carr's thesis [15], 


and they will be described there. 
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Assume 


THETAE 


Kitter=0 


POINTS 


FRONT 


Kitter=0 4 ITE 


Kitter=1 


89 


Compute 


Compute 


dy = 
2) to) 


90 


Compute 
yh) 


Compute 


x, (1) 


Compute 


(8 +8.) 


FRONT 


RADN2 


ok 


INDEX 


Compute 


y,(2) 1X, (1) 


Return 


POINTS 


Compute 


x, (J) 


Compute 


Y, (J) 


Compute 


Or (OO+6 4) 


Compute 


*, (9) +4 (9) 


Compute 


x, (J) 


a3 


FRONT 


Compute 


C setun _) 
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Compute 


x0,y0,z0 


N =-sin(a) 
x: 


N =cos(0) 
cos (a) 
N =sin(o) 


N_=-cos (a) 
Zz 


Change sign 


in N_,N 
5 gaat 4 
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INDEX 


Compute 


8,0,6' 


Compute 


KY pee Me 


Compute 


Compute 


(6546 .) WY gr 


Iloop=1 
5 


€=10 


RADN2 


6 j= (8, +6,) - 


Q 


Xn COS (89) 


RADI 


NO 


Compute 


8 (NEW) 
° 


Diff= 


I6 . (NEW) 
-8 (OLD) | 


ITloop>100 =ES 


NO 


Iloop=Iloop 
+1 


Compute 


peers peas 1S 
rx" Sry es 


Compute 


iP ,NP 
x 


fy’ 


INDEX 


100 


Compute 


p 


xr 


Compute 


Compute 


0,6,761 


ICHECK=50 


Compute 


102 


103 


Compute 


r (J) -geom. 


€ , (NEW) = 
8 


- 2 
OLD) ag 


zi 


NO 
Return 
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Iloop>90 


NO 


Iloop=Iloop+1l 


POINTS 


FRONT 


jsign=+1 


Compute 
Tei 


106 


108 


Compute 


' | ' 
a ee M 


109 


Compute 


K* Ty 
bai oe ne 


M° 
r 


Return 


Compute 


Kt sey M! 
ig If 1G 


P10 


Po 


7. COMPUTER PROGRAMS USER'S GUIDE 


Input Data 


BO. U5 Tae 


TR yen (15 
Nie Now No 


O, Grid, Ellnum 


Xar Div 


Shape, Squray, Ellips 


Options 
I = 4: R £ i i 
andes ear surface is a fixed cone 
(Capt. Carr's thesis, see ref. 
I =2 i i 
pace Front surface is a fixed cone 


(This report). 


Execution Commands 


: ae gee I Opening Commands 


LAXXXX (Write account no. instead of XXXX). 


Password 


GLOBAL TXTLIB FORTMOD2 MOD2EEH 


Fe ee Compilation 
FORTGI LENS (or GRIN) 


7.3.3 Run on Terminal cworke ‘2? 


FILEDEF 05 DISK DATA L 


FILEDEF 06 DISK OUT D (RECFM F BLOCK 132 PERM 


EXEC RUN LENS (or GRIN) 


XEDET OUT D 


se Ba 


(15]) 


7.3.4 HardsCopy 
FILE 
PRINT OUT D 
iwSeDe) GEL la PLOr 
LENSCOM PRINTER 
7.3.6 Change in Plot Routine 
XEDIT LENSPLOT FORTRAN 
7.3.7 Transfer Files to Another User 
(a) The Transferer 


CP SPOOL PUNCH CONT TO xxxxp {2) 


PUNCH LENS FORTRaN (3) 


(PUNCH GRIN FORTRAN (3) ) 


(3) ) 


CP SPOOL PUNCH CLOSE NOCONT 


(PUNCH DATA L 


CP SPOOL PUNCH OFF 
(b) The Receiver 
READCARDA *A * 
(c) Notes 
(1) WORKG is the name of the file of the 
commands in this section. 
(2) Account No. of the receiver. 


(3) Optional. 


Li 


8. COMPUTER PROGRAMS LIST OF SYMBOLS AND SUBROUTINES 


a ce a ce a a a ES 


8.1 Subroutines 


DATTA 


SHAPET 


SHAPEC 


INIT 


POINTS 


FRONT 


REAR 


CALTC 


NOSE 


ZOYO 


RADI 


SKEWT 


SKEWC 


INTENS 


IMAGE 


SUMM 


PLOTT 


INDEX 


RADN 2 


ITER 


DIRECT 


Initial data 


Computes shape of lens in Lae 


Computes shape of lens in I 1 


case 


Computes behavior of first ray. 

Computes coordinates of rear and front surfaces. 
Computes T,(9), I,'(J),¥,E at front surface. 
Computes n(J) ,I5(J) ,T5" (J) at rear surface. 
Computes thickness at each point for SHAPEC. 


Computes extrapolation of last points in I Ate 


case 


or coordinates shift in I 2. 


case 


Computes 2,'(IZ),y)'(1IZ) for skew rays calculation. 
Computes intersection of the ray with the designed 
surface. 


Skew rays calculation in I 


case 


Skew rays calculation in I 


case 


Calculates transmitted intensity of the rays. 
Calculates image plane spot diagram. 

Calculates spot diagram energy density. 

Save output for plot routine. 

Computes index of refraction in the lens (N5). 
Computes radius to rear (or front) surface from 
known radius to front (or rear) surface and 6_. 


e 


Iteration routine to find correct C. in SHAPET. 


Computes direction cosines of the ray at the rear surface. 
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Sue 


Se2el 


Se2ee 
x, (J) yy (9) 
Xo (J) +o (J) 
red) 
r (dg) 
XO, YO, ZO 
Kali y Noli dine 


Yom? “tm 


Program Symbols 


Data Symbols 
Focal length. 


The fixed surface cone half angle. 

The incident ray offset angle. 

Number of iterations. 

Maximum radius of fixed surface cone. 
Thickness of lens at the edge. 

Indices of refraction. 

inclination, of object plane. 

Size of grid increment. 

Number of ellipses. 

Z - coordinate increment for ellipse plot. 
Location of center symmetry in GRIN. 


The maximum change in the index of refraction 


(N5) 

Geometrical Symbo is 
Coordinates on front surface. 
Coordinates on rear surface. 


Radius to front surface. 

Radius to rear surface - 

Intercept of the ray with the front surface - 
Intercept of the ray with the rear surface - 
Intersection with the image plane. 

Direction Cosines Symbols 


Direction cosines of the front surface normal - 
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eee yi’ zi 
RRR, 
roe 7M’ 
Tex! fy!" £2 
NP eg NP ey NP ey 
Trex’! Try’ Trz 
Ved Pe a pa" 
8.2.4 


Direction cosines of 
Direction of coSines 


Direction cosines of 


ray. 


Direction cosines of 
Direction cosines of 
the ray at front surface. 
Direction cosines of the 


Direction cosines of 


refracted ray. 


Angle Symbols 


the rear surface normal. 
of the ray. 


the intemal refracted 


ErOnt racdLus. 


normal to plane of 


rear radius. 


rear surface external 


The angle symbols are defined in ficures 3 and 8. 
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9. PROGRAM LENS - LISTING 
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11. PRELIMINARY RESULTS 


Table 1. Comparison of RMS Spot size for a homogeneous lens 


with the RMS Spot size for a GRIN lens. 


+0.05 -0°.05 


SOT 
0.0057336 0.0069021 0.005286 


0. 

Get 0.027424 0.027646 0.027327 

OeZ 0.06277 0.06288 0.06304 

6.3 0.113431 O... 11216 0.113245 
| 0.4 0}. 18265 0.18287 0. 182 Fe, 
\ 


Note: DIV is the fractional variation of refractive index. 


a. is the angle of incoming parallel rays relative to lens 


12) 
axis; a> has dimensions of radians. The center of symmetry 
Xa was at -0.07 measured from the apex of the outer surface. 
The minus sign indicates that Xa is external to lens. 
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RESULTS FROM LENS 


(TCASE =.) 


A = 0.0, 0.2, 0.4 
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FOCAL LENGTH = 4.0000000 


R = 1.0000000 T = 


ALFAP = 0.0 


RAYS = 328 COUNT 


RMSRAD = 0.0000052 


ALPHA = 0.7853979 U = 0.0 1 
0-0500000 INDICES OF REFRACTION ARE? NL = 
SRID = 0.1000000 STATNA = 0.44171 GAMMA 
= 328 YCENTR = 0.0 SIGMAY = 0.0 


+ 


= 1008 


1.0000090 N2 


= 0.7850314 


= 1.5000000 


xC=0.0 


SIGMAZ = 0.0 


N3 = 1.0000000 


DIV=0.0 
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12. SUMMARY OF THE REPORT 

Lens design and ray tracing to determine optical 
aberrations requires efficient computer codes. A clearly 
structured computer code with efficient algorithms has been 
developed to design a lens incorporating gradient refractive 
index, GRIN. One surface of the lens, either the inside or 
outside surface, is a circular cone. 

The computer code includes lens design and calculation 
of optical aberration for a lens using a homogeneous material, 
i.e., for n equal to a constant. The homogeneous case is 
included for several reasons. The simpler case with n equal 
to a constant provides insight to the more complex case of 
variable n. Extensive computer codes require a check on the 
validity of the output. The homogeneous case provides a check. 
When the factor B is equal to zero in equation [2.9], the 
results for the two computer programs with and without GRIN 
should be identical. Needless to say, the results are identical 
for the programs reported here. 

The case of a lens with GRIN and with a circular cone for 
the outside surface has been programmed. The companion case 
with the circular cone on the inside is being developed by 
Carr [15]. Carr's program can be readily incorporated into 
this computer code. 

The program has been exercised for the case of 5% 


variation in refractive index for both negative and positive 


283 


values for B. The aim of the preliminary calculations 

was to verify operation of the computer program. More extensive 
calculations are required in order to understand the influence 
of various variables. As more experience is gained, the 
program should be combined with an optimization code to 

permit thorough and efficient lens design. 

The selection of a 5% variation in refractive index 
was based on current state-of-the-art in the production of 
GRIN materials. A wider variation in refractive index may 
be needed. 

The program is restricted currently to values of Xo 
satisfying the constraint xX, <¢ 9. Additional work is needed to 
permit designs using positive values of Xu: When Xo is 
positive, the angle & at some point along the ray interior to 
the lens has a critical value at which dr/dé@ changes sign. 


The critical value causes computational difficulties. 
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